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We have studied numerically the spectral correlations in 
a metallic phase and at the metal-insulator transition. We 
have calculated directly the two-point correlation function of 
the density of states R(s,s'). In the metallic phase, it is well 
described by the Random Matrix Theory (RMT). For the 
first time, we also find numerically the diffusive corrections 
for the number variance (Sn 2 (s)} predicted by Al'tshuler and 
Shklovskii. At the transition, at small energy scales, R(s — s') 
starts linearly, with a slope larger than in a metal. At large 
separations \s — s'\ 3> 1, it is found to decrease as a power 

law R(s,s') r c/\s - s'| 2 ~ 7 with c ~ 0.045 and 7 ~ 0.7, 

in good agreement with recent microscopic predictions. At 
the transition, we have also calculated the form factor K{i), 
Fourier transform of R(s-s'). At large s, the number variance 
contains two terms (Sn 2 (s)) = B (n) 1 +2tt K (0) (n) where K(0) 
is the limit of the form factor for t-»0. 



I. INTRODUCTION 

It is well established that the statistical properties of 
energy levels in an isolated metal present universal fea- 
tures characteristic of the chaotic systems QJ^]. At low 
energy scales, in the so-called ergodic regime, the correla- 
tions functions are well described by the Random Matrix 
Theory (RMT) ||. In the diffusive regime, i.e. when the 
electronic mean free path l e is smaller than the typical 
size L of the system, deviations from RMT occur Q when 
the energy scales become larger than the Thouless energy 
E c = T1D/L 2 . This energy is the inverse of the charac- 
teristic time to = L 2 /D for a particle to diffuse through 
the sample. D is the diffusion coefficient. On the other 
hand, in the localized regime, the correlations between 
levels are weaker and in the limit of an infinite system 
the statistics of energy levels becomes Poissonnian. 

It has been argued that the statistics of energy levels 
at the metal-insulator transition (MIT) is characterized 
by a third universal distribution which is a hybrid be- 
tween the Wigner and the Poisson distributions. This 
was clearly shown by a numerical study of the scaling of 
the nearest neighbor level spacing distribution P(s) ||. 
Several other numerical works have confirmed this idea 



f|s| — |l3|| - Moreover, by using diagrammatic calculations, 
semiclassical description and scaling arguments, it was 
found that long range correlations between levels exhibit 
a new power-law behavior whose exponent is related to 
the exponent v of the localization length |L4] [nj . 

There are several quantities which measure the fluctu- 
ations of energy levels. In the RMT, they depend only 
on the symmetry of the hamiltonian. If it is invariant 
under time reversal symmetry, the fluctuations are de- 
scribed by the Gaussian Orthogonal Ensemble (GOE) of 
random matrices ( /? = 1). When time reversal symme- 
try is broken, the spectrum becomes more rigid (Gaussian 
Unitary Ensemble, GUE, (3 = 2). These quantities are 
the following: 

• The number variance Y. 2 (E) 

Y> 2 {E) = {5n 2 {E)) = {n 2 {E))-{n{E)) 2 . (1) 

It measures the fluctuation of the number of levels 
n(E) in a strip of width E. The average (...) can 
be taken cither over different regions of the spec- 
trum or over several realizations of disorder. In the 
RMT, S 2 (i?) increases logarithmically with E. For 
E > A, it varies as T, 2 (E) ~ (2//3tt 2 ) hx(n) where 
(n) = E/A, A being the average interlevel spacing. 

• The distribution P(s) = (5(e — e$ + e»+i)) of the 
spacing s — e/A between consecutive levels. In the 
RMT, it is well described by the Wigner-surmise: 
P(s) oc s 13 exp(— cps 2 ), while, when there is no cor- 
relation between levels, it has a Poisson behavior: 
P(s) — exp(— s). 

• The two-point correlation function of the Density 
of States (DOS)Q: 



R(s,s') = \(p(e)p(e')) 
Po 



1, 



A < 2 » 



where po is the average DOS. Using the definition 
of the DOS, R(s) can be simply rewritten as: 



R(s = e/A) = (S(e - (e< - 
= 5(a) - 1 + 



e0)>"l 
P(n, s) 



E 

n 



(3) 
(4) 



1 This function is usually called K(e) However, we have 
chosen to follow here the notations of refs. [I14l-|l9i 
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where P(n, s) is the distribution of distances s n be- 
tween levels separated by n other levels By 
definition P(s) — P(0, s). The number variance 
can obviously be written in terms of this two-point 
correlation function: 



Y?{E) = / R(s - s')dsds' (5) 

Jo Jo 

(n) 

(6) 



2 1 ((n) -s)R{s)ds. 
la 



The form factor K(t), Fourier transform of R(s): 



K(t) = — I R(s) exp(ist)ds . 

Z7T 



(7) 



The interest of this quantity is that is can be di- 
rectly related to some characteristics of the classical 
motion of the diffusive particles J22|. 

Although the number variance is most frequently used 
in the literature, it is not always the most appropriate 
quantity to describe the correlations, since it is a dou- 
ble integral of the DOS-DOS correlation function R(s). 
Thus the behavior of this quantity at an energy scale E 
depends on the behavior of R(s — e/A) for all energies e 
smaller than E. That is why in this paper we have cho- 
sen to study R(s) as well as its Fourier transform K(t) 
directly. Although analytical predictions have been given 
for the asymptotic behavior of these quantities (large en- 
ergies or small times) at the MIT ]T^-pj[|, nothing was 
known yet about the short range behavior of the correla- 
tions. It is one of the goals of this paper to study these 
correlations. In section 2 we recall what is their behav- 
ior in the metallic regime and compare with numerical 
experiments. In section 3, we analyze the two point cor- 
relation function at the MIT. Section 4 is devoted to a 
discussion of level spacing distributions, and finally, in 
section 5, we present our conclusions. 



II. SPECTRAL CORRELATIONS IN THE 
METALLIC REGIME 

Efetov has shown that at energy scales smaller than 
the Thouless energy E Cl and in the limit where the di- 
mensionless conductance g — E c /A is much larger than 
1, i.e. far from a MIT, the two-point correlation function 
R{s) is given by its expression in the Gaussian RMT 
[jlj . This was confirmed numerically in the framework of 
the tight-binding Anderson model with diagonal disor- 
der: the number variance and the distribution P(s) were 
found to be very well described by the Wigner-Dyson 
predictions 

We first show numerical results for the two-point cor- 
relation R(s), fig.|l], and its Fourier transform K(t) in the 



metallic regime, fig.|[ It is seen that the RMT result is 
a very good description of these correlations. They were 
obtained by exact diagonalization of a tight binding An- 
derson hamiltonian with on-diagonal disorder via a Lanc- 
zos routine. The diagonal elements are box-distributed 
around zero with a width W, in units of the transfer in- 
tegral t. We only use eigenvalues which lie in a central 
energy strip whose width is chosen such that it contains 
approximately half the number of all eigenvalues. Even 
though the average density of states (DOS) is almost con- 
stant in this region, special care has to be taken in un- 
folding the spectrum. Then we calculate R(s) defined 
in eq.[|. Thus, i?(e) is the distribution of differences of 
eigenvalues in the energy interval AE used for averaging. 
In addition, we average on different disorder configura- 
tions. The self-correlation of the levels is not included 
in the numerical calculation of R(s). A 5 function has 
to be added to the result. Using rather small systems 
but many disorder realizations (8x8x8 sites and 1000 
disorder realizations for e < 5), we are able to obtain a 
large amount of eigenvalues and correspondingly rather 
smooth curves for i?(e). As shown in fig.[l]b, the numeri- 
cal precision is even high enough to allow for the obser- 
vation of the tiny oscillations in R(e) in the GUE case 
where R(s) = S(s) — sin 2 (7rs)/(7rs) 2 . 

Let us briefly recall some important features of these 
correlation functions. It is helpful to use a semiclassical 
argument developped by Argaman et al. ]22| | in which the 
function K(t) is directly related to the return probability 
Pit) after a time t for a classical diffusive particle. More 
precisely: 



K(t) = 



2tP(t) 
4tt 2 /3 



(8) 



The factor 2 in the numerator arises from interference 
effects between time reversed paths, in the absence of 
magnetic field or flux. This relation is valid for classical 
times t much smaller than the Heisenberg time th — fi/ A 
For time scales larger than td = Ti/E c , the diffusion 
is homogeneous so that P(t) is constant. This gives a 
linear behavior for Kit) as observed in fig.0. For t ^> 
th, K(t) has to saturate. In the energy space, K{e) 
has first a 8 peak at the origin which describes the self 
correlation of the levels. Fourier transform of this peak 
implies that Kit) — > 1/2tt for large t. For small energies, 
there is a relation between the function i?(s) and the 
distribution of interlevel spacings P{s): R(s) = S(s) — 
1 + Pis) so that indeed i?(s) starts linearly at small s in 
the absence of a magnetic field. At larger energies, -R(s) 
varies like — l/(/?7r 2 s 2 ). Unless otherwise specified, we 
will now consider the case (3=1, i.e. no magnetic field. 



2 For a discussion on the validity of eq.^ and characteristic 
times, see ref. [E21 
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For time scales smaller than td, the diffusion is non 
homogeneous and the classical return probability varies 
as 



P(t) 



V 



(A-KDty/ 2 



(9) 



so that K(t) scales like i 1 ^/ 2 . 
transform R(s) is given by: 



For e ^> E Cl its Fourier 



R(s) cx 



1 



s \ d/2 
\ 9 ) C ° S{ - ] 



(10) 



and the number variance behaves like {E/ E c ) d / 2 for E > 

This " Al'tshuler-Shklovskii" power-law behavior had 
not been seen numerically yet . This is because there 
is a long cross-over regime around the energy scale E c , 
and the power law is expected to be visible for very large 
energy scales E E c . In the intermediate regime E > 
E c , the quantization of the diffusion modes has to be 
taken into account: 



P(t) = ^exp(-iVi), 



(11) 



Using a large system with small E c , we have calculated 
£ 2 (£) numerically and found the onset of the power law. 
S 2 (£') starts to be proportional to {E/E c ) 3 / 2 only for 
energies E > 100£' c (see fig.|j). If one takes into account 
the fact that for a finite system of size L 3 the condi- 
tion n <§; L 3 should be fullfilled in order to neglect finite 
size effects, and if one wants to observe the power law 
over at least one decade of energy, the size of the system 
should be at least of the order L 3 ~ 10 4 . On the other 
hand, the discrete sum describes the numerically calcu- 
lated Y?(E) very well, and a fit to this sum with E c as 
only fitting parameter allows for a precise determination 
of the Thouless energy of the sample (fig-|3|) ■ 

It turns out that this contribution of the finite-q diffu- 
sion modes is hardly visible on the two-point correlation 
function K{s). For a metal in three dimensions, K(s) 
varies at large energies positively like +l/^/s instead of 
— 1/s 2 for RMT, meaning that the levels tend to attract 
each other at large energies p3fl . However, these devia- 
tions from RMT are very small, less than a few percent. 
We have not been able to see them on the numerically 
calculated K(s), see figjl]. For the same reason, the devi- 
ations from the linear behavior of K(t) for t < t& are not 
visible. Moreover, in the metallic regime that we have 
studied (g 3> g c ), we did not see the finite conductance 
deviations at small energy scales s ~ 1 predicted recently 



III. THE SPECTRAL TWO POINT 
CORRELATION FUNCTION AT THE 
TRANSITION 



We now turn to the description of the correlations at 
the MIT. The first theoretical predictions of spectral cor- 
relations at the MIT different from those in the metallic 
regime were given by Al'tshuler et al. §. They pre- 
dicted a purely linear behavior for the spectral rigidity, 
Y?{E) = A(n) for {A) » 1, where the constant A was 
estimated to be of the order of 0.25, and thus consider- 
ably different from the behavior in the insulating regime 
where A is unity. Later, Shklovskii et al. argued that 
the level statistics at the MIT should be completely uni- 
versal (|]. They predicted a level spacing distribution 
which should be a hybrid between a Wigner-Dyson type 
behavior for small s, that is P(s) cx s for s — > 0, and 
a Poisson type for s 3> 1: P(s) oc exp(— as). Then 
Kravtsov et al. brought new dynamics into the field by 
showing that the two-point correlation function R(s) has 
a new long range power-law behavior: instead of vary- 
ing as — 1/((3tt 2 s 2 ) like in the Wigner-Dyson regime, it is 
predicted to obey (wl : 



R(s) 



Cd/3 
„2- 7 



(12) 



The non-universal constant Cdp depends on the dimension 
and on the symmetry of the ensemble, and 7 is given by 
7=1 — l/(yd), where v is the critical exponent of the 
localization length £. If one uses the value v ~ 1.3—1.5 at 
the MIT p6fl , the exponent 7 is of the order of 7 ~ 0.75. 

This new power-law behavior can be derived semiclas- 
sically from K(t) along the same lines as in the metallic 
regime, but now with anomalous diffusion p^l7|. Near 
the MIT, the dimensionless conductance g = E c / A scales 
with distance A like g(A) = g c [l + (A/^) 1 /"], where v is 
the exponent of the localization length £(g). Since A 
scales like A~ d and E c like A -2 , this implies that the dif- 
fusion coefficient has to be size dependent. For A £ it 
has to vary as 



D(A) cx ^A 2 ~ d 
Po 



(13) 



To this length dependent D(A), one can associate a time 
dependent diffusion coefficient D(t), 



Din v (^) J "'/ ' ; 



A I q r t N 
1 + - — 



(14) 



so that according to eqs. || and K{t) is constant up to 
a time dependent correction: 



1 - (^(ScAi) 1 ^ 



(15) 
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The semiclassical expression (||) is known to be equivalent 
to the two-diffuson diagrammatic result of Al'tshuler- 
Shklovskii 0, now with an energy dependent diffusion 
coefficient. Aronov et al. |ll| argue that at the tran- 
sition, there are many more diagramms with the same 
structure so that the correct prefactor of the power law 
is not known, even in sign. This structure is supposed to 
be correct for t <C th- By Fourier transform, one has the 
tail of R(s) 



R(s) oc - 



l/u 



l + l/vd 



(16) 



At the transition £ oc L so that one finds the size in- 
dependent power law for the tail of the two-point corre- 
lation function, eq.^. The integration of the power-law 
tail gives a contribution to the number variance S 2 (_B) 
varying as (n) 7 at large energies. The existence of a lin- 
ear term in £ 2 was first ruled out because of a sum rule 
jl4| , [l5f — a statement later discussed and contested by 
Aronov and Mirlin (T^JUl- 



Finally, Sears and Shore |10| introduced the function 
r(z) = S 2 /(n), which they computed numerically. With 
the hypothesis that the scaling behavior of Y(z) should 
only depend on the combination of parameters z = 
(n)(£/L) d , they were able to reproduce the term (n) 1 in 



E 2 with the same 7 as introduced in |14|. However, they 
indeed found an additional linear term. Their numerical 
results are in agreement with T(z) = 0.30 + 0.22z~ 2 / 9 for 
a rather large range of disorder and for system sizes up 
to 16 x 16 x 16. 

All in all, a consensus now emerges that for (n) 3> 1, 
the number variance should contain two terms 



E 2 (E) = A(n)+B( n y. 



(17) 



Although the power law originates directly from the 
long range behavior of the two-point correlation function, 
the origin of the linear term is more subtle. It appears 
as a constant of integration which is a global measure of 
the interaction between eigenvalues integrated over en- 
ergy. More precisely, since by definition (eqj^), R(s) is 
the second derivative of S 2 (s), the coefficient A of the lin- 
ear term is simply given by A = J_ oo R(s)ds = 2irK(0). 
The absence of a linear term in T, 2 (E) in the RMT im- 
plies that the "sum rule" R{s)ds = is fulfilled. 
The coefficient A contains informations about all energy 
scales, in particular low energy scales. However, until 
now, most studies concentrated on the long range be- 
havior of (5n 2 ). In the present work, we focus on the 
direct calculation of the two-point correlation function 
and of the form factor both for s < 1 and s ^> 1. Using 
the method described above, we have calculated R(s) for 
rather small system sizes but using many disorder real- 
izations for W = W c = 16.5. In the energy range we 



considered, finite size effects can be neglected, which we 
checked by using different sizes. In figj^ we show the uni- 
versal results obtained from systems with different sizes 
in an energy range < E < 2. In this small energy range 
we used system sizes uptol6xl6xl6 sites and up to 
1000 disorder realizations. Obviously, the correlations 
are weaker than in the GOE and qualitatively in agree- 
ment with the analytical result of Kravtsov and Mirlin 
p4j which predicts a correction to the small energy slope 
increasing like l/<? 2 when g decreases. 

The linear variation of R(s) at small s is consistent 
with the behavior of P(s) at small s (fig.||). More pre- 
cisely R(s) = 6(s) — 1 + P(s). This means that all the 
P(n, s) for n > 1 start at least quadratically for s — > 0, 
like in the RMT §J. 

For practical purpose , we have fitted R(s) for s < 5 
with 



R(s) = 8(a) - a * ex P(- s M) (18) 
i=i 

The found coefficients and Sj are shown in table |. 



i 


1 


2 


3 


4 




0.74 


0.23 


0.02 


0.01 


Si 


0.17 


0.51 


2.85 


3.73 



TABLE I. Numerical coefficients that describe the short 
range correlations (see Eq. (|l8|)). 



It is seen in fig.|| that for e > 2, R(e) is very close to 
zero and the correlations become very difficult to find. 
In order to detect correlations at larger energies, we have 
diagonalized a total of 132 000 systems, the system sizes 
ranging between 6x6x6 and 8x8x8 sites (see table 
0). 



System 


d 


c 


2-7 


6x6x6 


30 000 


-1.17 


0.047 


6x7x8 


45 000 


-1.12 


0.039 


7x7x7 


20 000 


-1.32 


0.046 


7x8x9 


19 000 


-1.08 


0.035 


8x8x8 


18 000 


-1.14 


0.037 


average 


total: 132 000 


-1.17 ±0.08 


0.041 ± 0.005 



TABLE II. Numerical values for the power-law behavior 
R(s) oc — c/s 2-7 for five different geometries, d is the number 
of disorder realizations. 
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Analyzing the central half of the spectra (totaling 
about 23 million levies), we have been able to extract the 
behavior of R(s) for s < 25. In fig. |] we show the average 
correlation function obtained from all systems together 
with the best fit of the tail to a power law, given by 



R(s) 



0.041 

„1.17 



(19) 



in good agreement with the predicted eq.[l2| The compar- 
ison of the fitted variables for the five different geometries 
used allows to estimate the remaining statistical error 
for the exponent to 1.17 ± 0.08 and for the prefactor to 
0.041 ± 0.005 (simple standard deviations, see tabic ||). 
Rather than by statistical errors, the calculation of R(s) 
for even higher energies is limited by the precision of the 
unfolding of the spectra: We were able to obtain (/9(e)) 
constant with a precision of about 10~ 4 . For s ~ 30, R(s) 
becomes comparable with this value and thus starts to 
fluctuate around zero. This explains the sharp drop in 
fig](] at these energies and might cause a further system- 
atic error for s smaller but close to 25. We estimate the 
total error of the exponent to be less than about 20%. 

We now turn to the discussion of the sum rule. As we 
have seen above, the coefficient of the linear term in S 2 
is the integral of R(s). It can be written as 



A = l 



R(s)ds 



(20) 



o+ 



In the case of the Wigner-Dyson statistics, the weight 
/ 0+ R(s)ds — —0.5, exactly compensates the weight 1 of 
the 5 function at the origin. In the case of a Poisson 
statistics, R(s) is zero for s ^ 0, so that the S function is 
not compensated and there is a linear term with slope 1 
in £ 2 . The situation at the MIT is intermediate. 

We see here that the prefactor is modified by the short 
range behavior of R(s). Inserting the numerical values 
from the fit, we find that A(s) = 1 — 2 f Q+ R(s')ds' is 
about 0.37 for s — 5 and 0.22 for s — 25. Long range cor- 
relations reduce A further: assuming that one can extend 
the obtained power law for R(s) up to much higher ener- 
gies, one can even estimate A(s) for very large energies. 
It gives A ~ 0.17 for s = 100, A ~ .09 for s = 1000 and 
... A ~ for s — > oo. The sum rule seems to be obeyed 
for energies up to infinity. This result constradicts the 
results of refs. || and Jl0| who predicted A(oo) ~ 0.25 
and A(oo) ~ 0.30, respectively. 

We finish this study of the two point correlation func- 
tion by showing in fig.0 the numerically calculated stru- 
ture factor K (t) as a function of time, including ener- 
gies up to 25 A. As expected from the above semiclassi- 
cal argument, it tends to a finite value for small t (but 
still t has to be larger than the collision time). This fi- 
nite value is related to the slope of Yj 2 (E) for E — > oo: 
K(Q) = limB^oo^E 2 \E)/dE. In our plot, K(t) seems 
to converge to the value K (0) ~ 0.4. However, in the 



numerical evaluation of K{t) one has to use a cut-off in 
energy. Using E ~ 25, our data for K are valid only 
down to t ~ 0.04, so that we are not able to calculate 
properly K(Q). For the moment, we can only conclude 
that a time dependence of the form 



k (t) = k(o) + m 1 -' 1 



(21) 



is compatible with our data (b > 0) . 

For the same reason as in the metallic case (self corre- 
lation of levels), k{t) has to saturate at k(t) — 1/2-k for 
large t. 



IV. LEVEL SPACING DISTRIBUTIONS AT THE 

MIT 

For completeness, we briefly comment on the results 
concerning P(s). Shklovskh et al. || argued that at 
large s, it must follow a "modified" Poisson distribu- 
tion of the form exp(— as) while Aronov et al. predict 
a tail of the form exp(— Ad(3s 2 ~ 1 ). This tail reflects the 
behavior of the two-point correlation function at large s. 
This is found by describing the thermodynamics of a one- 
dimensional gas with a pairwise interaction consistent 
with the above two-point correlation function |ll|[l8|,|l9| . 
There is no complete consensus yet on the numerical de- 
scription of P(s) at the MIT HQ but it seems that it 
can be globally well fitted by an expression of the form 
i?i(7)s exp(— i?2(7)s 2-7 ), the constants B\ and E>2 being 
fixed by normalization of P(s) pi]Jl3[ |. Quite recently 
Kravtsov and Lerner have constructed a plasma model 
with the appropriate power law interaction in order to 
get the correct dependence of the two-point correlation 
function (|l^). They have been able to relate the coeffi- 
cient £?2 to the prefactor c of the power law . They find a 
value of c ~ 0.03 in reasonable agreement with our result. 

One may wonder if the knowledge of the probability 
distribution of distances between non-neighboring levels 
could add more useful information on the long range cor- 
relations. For this purpose, we have studied the function 
P(n, s) defined above. These functions are known for the 
RMT regime Q|. In the Poisson limit, they simply vary 
as: 



P(n,s) 



s'"e 



(22) 



It may be much more useful to use these quantities be- 
cause, for large s, R(s) and P(s) are very small and diffi- 
cult to calculate. On the other hand, for large distance s, 
P(n, s) is very different from the RMT or Poisson cases. 
Fig.|| shows the successive P(n, s) at the MIT, and they 
are compared to the corresponding functions in the Pois- 
son and RMT cases. For s ^> 1, P((s),s) characterizes 
much better the differences between the statistics at the 
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MIT and in the metal than P(0,s). Unfortunately, we 
do not know any prediction for the behavior of P(n, s) 
for large s, n at the MIT. 

V. CONCLUSIONS 

We have presented results of the direct numerical cal- 
culation of the spectral two point correlation function 
in metals and at the metal-insulator transition. In the 
metallic regime the results are in very good agreement 
with the RMT predictions. We have been able to find 
numerically the Al'tshuler and Shklovskh regime where 
the diffusion is power-like. At the metal-insulator tran- 
sition, the short range behavior of the two point corre- 
lation function is clearly modified. The short range part 
of the correlation is weaker and is followed by a power 
law behavior with an exponent which is consistent with 
the picture of anomalous diffusion. At the transition, 
R(s) is size independent. This clearly refutes a predic- 
tion |l^] that the power-law contribution should decrease 
as L — > oo. These short range correlations give rise to a 
non-trivial linear term in the spectral rigidity Y?(E) in 
addition to a power-law term which is the signature of 
the long range correlations. By integration of R(s), we 
conclude that the prefactor of this linear term may be 
much smaller than previously anticipated and may even 
be zero. 

Recently, different models of interpolating ensembles of 
matrices have been proposed p4]-|37|] . Our result puts a 
constraint on the possible choice of ensembles since they 
should produce a form factor which varies at small times 
like eq.pl|. 
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Kravtsov, I. Lerner, M.L. Mehta and B. Shklovskh. Nu- 
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puter of IDRIS (Institut du Developpement et des 
Ressources en Informatique Scientifique, Orsay). 



FIG. 1. Correlation function R(s) vs s in the metallic 
regime (8x8x8 sites, W = 4) a) with no external flux, b) for 
4> = 0.2500. The continuous lines are the predictions of RMT. 
The energy E is measured in units of A, the mean interlevel 
spacing. We have omitted the 5 function at the origin 

FIG. 2. The Fourier transform K(t) versus time t in the 
metallic regime a) with no external flux, b) for <f> — 0.25<^o. 
The continuous lines are the predictions of RMT. t is mea- 
sured in units of h/A. 

FIG. 3. Number variance E 2 (i5) in the metallic regime 
for a system with 20 x 20 x 20 sites and E c ~ 2.5A (empty 
circles). This figure shows the cross-over regime between the 
ergodic (RMT) regime (dotted line) and the E 3 ^ 2 regime (full 
circles). Deviations from the RMT behavior become visible 
for E of the order E c . The cross-over regime is well fitted 
by the discrete sum over modes eq.[ll] (full line). The onset 
of the asymptotic E 3 ^ 2 behavior is visible for our finite size 
system. For larger energies E 2 (_B) decays again due to finite 
size effects. 

FIG. 4. The short range behavior of the correlation func- 
tion R(s) at the metal insulator transition, obtained from 
systems with sizes 8x8x8 (dots), 14 x 15 x 16 (stars), and 
16 x 16 x 16 (diamonds). We used up to 1000 disorder re- 
alizations. The continuous line is the RMT result (metallic 
regime). We have omitted the 8 function at the origin 



FIG. 5. The correlation function R'(s) = R(s) + 1 - S(s) 
together with P(s). The low energy behavior is identical. 

FIG. 6. The long range behavior of the correlation function 
—R(s) at the metal insulator transition on a logarithmic plot 
(full circles), together with the best fit to a power law (full 
line). The function was obtained from 132 000 systems with 
sizes between 6x6x6 and 8x8x8 sites (see table [□]) 



FIG. 7. Form factor K(t) at the metal-insulator transi- 
tion. The continuous line shows the RMT result which fits 
the metallic behavior. 



FIG. 8. P(n,s) in the metal (empty circles), for Poisson 
(full lines), and at the transition (full circles). 
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